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ABSTRACT 

The nonlinear interaction between two oblique three-dimensional Tollmien-Schlichting 
(TS) waves and their induced streamwise-vortex flow is considered theoretically for an 
incompressible boundary layer. The same theory applies to the de-stabilization of an 
incident vortex motion by sub-harmonic TS waves, followed by interaction. The scales 
and flow structure involved are addressed for high Reynolds numbers. The nonlinear 
interaction is powerful, starting at quite low amplitudes with a triple-deck structure for the 
TS waves but a large-scale structure for the induced vortex, after which strong nonlinear 
amplification occurs. This includes nonparallel-flow effects. The nonlinear interaction is 
governed by a partial-differential system for the vortex flow coupled with an ordinary- 
differential one for the TS pressure. The solution properties found sometimes produce 
a break-up within a finite distance and sometimes further downstream, depending on 
the input amplitudes upstream and on the wave angles, and that then leads on to the 
second stages of interaction associated with higher amplitudes, the main second stages 
giving either long-scale phenomena significantly affected by nonparallelism or shorter quasi- 
parallel ones governed by the full nonlinear triple-deck response. Qualitative comparisons 
with experiments are noted. 
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1. INTRODUCTION 


The idea of significant interactions arising between oblique Tollmien-Schlichting (TS) 
waves and induced, or incident, streamwise vortices has been circulating for some con- 
siderable time, as a possible or even a probable early stage in transition to turbulence in 
boundary layers and duct flows. In the boundary-layer context there is some experimental 
evidence of such interactions occurring in practice (Aihara et al., 1965, 1969, 1981, 1985; 
Tani and Sakagami, 1962; Bippes and Gortler, 1972; and possibly Klebanoff, Tidstrom 
and Sargent, 1962), as well as supporting computational work (see Wray and Hussaini, 
1984 and Spalart and Yang, 1986, the latter showing the emergence of oblique waves from 
random initial disturbances), for flat or curved surfaces. The theoretical idea is basi- 
cally this in principle: if two low-amplitude TS waves represented by (ax,±/?x,Dx) are 
present, i.e. each with a streamwise wavenumber a.i and a frequency fix, but with oppo- 
site spanwise wavenumbers ±/?i, then nonlinear inertial effects produce the combination 
(«i — ai,/?i + — fix) at second order, among other contributions, i.e., the standing 

wave or streamwise vortex flow (0, 2/?i,0) is induced. Equally, the combination of the 
vortex (0, 2/?i,0) and one TS wave (ax, ft, Ox) provokes the other TS wave (ax, — /?i, fix) 
also. So nonlinear interaction takes place between the three-dimensional (3D) TS waves 
and the vortex. Such interplay may start from a single 3D TS wave dominating at first, 
upstream, with a minor amount of the other TS wave and hence of the induced vortex, 
or it may be initiated by an incident vortex upstream, with the TS waves then developing 
at first as small disturbances of the boundary-layer-plus- vortex motion. Both situations 
are observed experimentally or computationally. In either situation the ensuing nonlinear 
interaction downstream may be expected to have the same character to some extent, al- 
though possibly with considerably different end results, depending on the input conditions 
upstream. There is a possible connection here in fact with the use of vortex generators in 
laminar-flow control (Bushnell, 1984) in boundary layers. 

Little progress has been made, however, in actually converting the above notion into 
nonlinear TS/vortex interaction equations in a rational theoretical way consistent with the 
Navier-Stokes equations. The only such theoretical study appears to be that of Hall and 
Smith (1987, see also references therein), but that applies to channel-flow interactions. A 
major difficulty exists in moving on to the boundary-layer context, namely to make al- 
lowance consistently for the growth or decay of the basic steady boundary layer, say on an 
airfoil, i.e. to incorporate nonparallel-flow effects. These effects are known to be substan- 
tial in Gortler- vortex development on a curved surface (Hall, 1982, 1988 and references 
therein) and in TS waves to some extent (Smith, 1979a and references therein), and like- 
wise they are potentially significant in the long TS-induced vortices or incident vortices of 
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present concern for a flat or curved surface. A second major task is to pin down the main 
scales, and hence the flow structure, involved in a TS/vortex interaction in a boundary 
layer, and the scales are also affected by the nonparallelism present. Our aim here is to 
incorporate these two aspects in a self-consistent theory to begin tackling boundary-layer 
TS/vortex interaction. Mention should be made at this stage of the theoretical approxi- 
mations by Benney and Lin (1960), Antar and Collins (1975), Srivastava and Dallmann 
(1987), among others, all of which are of undoubted interest in their attempts to address 
various levels of weak interaction, but all of which are in the end inconsistent with the 
Navier-Stokes equations. This is due much to the nonparallel-flow effects acting at finite 
Reynolds numbers: see also later comments. 

The current approach takes advantage of the feature that the Reynolds numbers (Re) 
of interest in reality are very large. So Re is taken as a large parameter throughout. 
Previous discussions (e.g., Smith, 1979a, b; Hall, 1982; Hall and Smith, 1984; Ryzhov and 
Zhuk, 1986; Smith and Stewart, 1987) show that the large-i?e approach is in fact the only 
way to obtain a consistent theory for boundary-layer instabilities and transition, apart 
from large-scale computations of the 3D unsteady Navier-Stokes equations of course, and 
thus, in a sense, there in no competing theory as far as we are aware. Further, as a bonus, 
the comparisons with experiments and computations presented in the above papers, when 
combined with the numerical interpretations of Smith, Papageorgiou and Elliott (1984), 
tend to verify the practical value of the current type of approach at the Reynolds numbers 

S' 

of real concern. 

The scales and the flow structure involved in the TS/induced- vortex interaction in an 
incompressible boundary layer are examined in sections 2-4 below. They are based only 
loosely on our channel-flow study (1987) and have to accommodate the shortness of the 
triple-deck streamwise length scale in the typical 3DTS wave compared with the longer 
scale of the typical vortex affected by nonparallelism (Smith, 1979a, b; Hall and Smith, 
1984). Whether the nonparallelism matters in the long run or not then remains to be 
seen, although clearly it can continually destabilize the boundary layer and bring in the 
influence of external pressure gradients for instance. In any case, the scales reflect the 
property that one of the strongest vortex effects arises from the slow decay of the 3DTS 
waves’ amplitudes at the edge of the viscous lower deck or critical layer. Thus the spanwise 
TS velocity (it;) there decays as the inverse of the normal distance y from the surface, and 
so the inertial response at second order is like y~ 2 . This must balance the main viscous 
force oc d 2 w/dy 2 of the induced vortex, implying that the vortex w-component grows 
logarithmically with y. As a result, substantial vortex motion occurs well outside the 
viscous TS layer. That substantial motion spreads way out across the boundary layer and 
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beyond, although its dominant nonlinear interaction effects with the TS flow turn out to 
spread no further initially than a thin buffer zone lying just outside the lower deck but 
within the boundary layer: see Figure 1 later. The scales can be derived in a number of 
ways, we note, for example as in Section 2 or by starting from 3D triple-deck theory and 
working outwards, as it were. The resulting flow structure found captures the original 
notion of TS/ vortex interaction in the form 

E X E? = E 3 (1.1) 

where E lti represent the oblique TS waves and E 3 represents the induced vortex, specifically 
Ei i = exp[*(aX ± (3Z/2 — flT)] and E 3 = exp(i/3Z) in the later notation. The TS-squared 
forcing of the vortex motion takes place both in the buffer layer and in the lower deck 
whereas the vortex reaction back on the TS waves is concentrated in the lower-deck flow. 
The production of the spanwise subharmonics E\ t i here is in line with Aihara and Koyama’s 
(1981) experimental findings in vortex break-up. 

The interaction equations are summarized and addressed numerically and analytically 
in Section 5. Certain of the nonlinear properties found appear to agree qualitatively 
with the Srivastava and Dallmann (1987) approximations a little and with the Aihara 
and Koyama (1981) experiments, but others do not. A point here with regard to the 
experimental comparisons is that the theory, being a first step (see earlier comments), 
deals with a first nonlinear-interaction stage occurring at low amplitudes. Second stages 
corresponding to higher amplitudes can arise either as a downstream development from 
the present stage or as a consequence of increased input amplitude upstream. Second 
stages are indicated directly by the current analysis, as the further discussion in Section 6 
describes, and the second-stage interactions can have shorter or longer streamwise extent, 
depending on the input conditions and the spanwise wavenumbers. 

The velocities «<»(«, v,t«), corresponding Cartesian coordinates ioo{x f y,z) and pressure 
Poo^loP are used, based on the airfoil chord 4o, free-stream speed u in the x-direction 
and the incompressible fluid density p <», with time written as u^L^t. To fix matters, we 
may consider usually the flat plate (y = 0 for 0 < x < 1) with zero pressure gradient, 
giving a basic Blasius boundary layer with normalized skin-friction factor A = 0.332 Xq at 
the typical station (x,z) = (x 0 , z 0 ) under consideration. The extension to nonzero pressure 
gradients, for example, is referred to in Section 6, and the extensions to nonzero Gortler 
numbers for curved surfaces and to the compressible regime are also of much interest. 
The terms A*,, p a introduced later are higher-order corrections of the basic velocity profile. 
Also, Re = ttoo^ 0 «w( = £ ~ & ) i s the global Reynolds number, with being the kinematic 
viscosity of the fluid, the superscript * denotes a complex conjugate and r,i stand for the 
real and imaginary parts of a quantity. 
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2. THE SCALES AND THE MAIN-DECK RESPONSE 


There are several ways to arrive at the scales controlling the 3D TS/induced-vortex 
nonlinear interaction of present interest. One is as follows: see also Fig. 1. Linear 
or nonlinear 2D and 3D TS waves are viscous-inviscid interactive phenomena, governed 
primarily by the triple-deck structure and hence by the unsteady interactive boundary- 
layer equations holding in the lower deck which has y-extent 0(e 5 ) near the wall (Smith, 
1979a). If the 3DTS amplitudes are comparatively small, say of order h relative to the fully 
nonlinear sizes, then inertial effects force a vortex motion at relative order h 2 . The span wise 
inertial effects ( udw/dx , etc.) in particular decay only slowly, however, like Y~ 2 , with 
scaled distance Y from the wall, due to the algebraic decay of the 3DTS velocities into their 
asymptotes, and so the spanwise velocity ( w v ) of the induced vortex grows logarithmically 
like tnY (Hall and Smith, 1984) since the vortex response is predominantly viscous here. 
This growth is then damped down further away from the wall where Y ~ <5(» 1) say 
in a buffer zone and the relatively slow streamwise variation of the vortex takes effect, 
introducing the comparatively slow inertial operator Yd-^ [from the boundary-layer shear] 
to combat the viscous one d\. So is 0(<5 -3 ). There w v is of relative order h?ln6, 
provoking by continuity a streamwise velocity u v of relative order h 2 8 z in8 typically, in the 
vortex, which therefore alters the basic shear by a relative amount of the order h 2 8 2 Cn8. 
The relative influences of the slow streamwise modulation ~ 8~ 3 and the change in mean 
shear ~ h 2 8 2 tn8 above are in balance if 8 ~ h~%, to within a logarithmic factor, and 
this balance is central to the nonlinear evolution of the provoked vortex flow and the TS 
disturbance. Another feature that needs to be considered is the effect of nonparallelism in 
the basic flow, due to the boundary-layer divergence among other things. Nonparallelism 
affects the TS and vortex interaction over unsealed streamwise lengths typically of order 
6 ~ 3 since that is the characteristic amount by which the critical linear conditions are 
disturbed in the vortex-TS interaction above. On the other hand, the unsealed length 
associated with the streamwise modulation of the TS wave or vortex is 0(e 3 £ 3 ) because 
the triple-deck length scale is 0(e 3 ). Hence nonparallel-flow effects are accommodated in 
the balance e z 6 z ~ 8~ z , giving 8 ~ e*, which is large as assumed initially. Nonparallelism 
would be negligible for 8 <C e - *, but the stage to be addressed here, namely 

h ~ e*, (2.1 a) 

is possibly a crucial one in our opinion as it marks the intrusion of significant nonparallel- 
flow effects and hence of more global properties into the nonlinear TS-vortex interaction. 
(We remark here that the scales in (2.1a) can be inferred from Hall and Smith (1984), and 
also if the induced vortex is absent the disturbance size required to bring in nonparallelism 
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is raised to h ~ e*). Whether the nonparallelism thus captured by (2.1a) actually has 
a substantial impact on the interaction or not can only be decided by the solution prop- 
erties of the resulting nonlinear governing equations. One inner limit of these equations 
reproduces the parallel-flow version of the TS-vortex interaction in any case. With the 
suggested scales (2.1a) holding, then, the scene is set for a multi-scaling in which 

— ► e d x — * e 3 dx + £ (2.l!>) 

where (t,x — xo, z — zo) = (e 2 T,e 3 X,e 3 Z) define the main 3DTS time and length scales, 
x — x 0 = erf X defines the slower modulation scale of the TS waves and the induced- vortex 
flow, and four y-scales come into operation: 

y = e 3 y', e*y, e*y, e 5 Y. (2.1c) 

The e 3 ,e*,e 5 levels are associated with the three decks of the triple-deck TS structure, 
while the extra level e* is necessary for the determination of the forced vortex motion 
because of (2.1a). The scales (2.1a) also indicate that the amplitudes involved in nonlinear 
TS-vortex interaction are quite small, initially at least, with a typical TS pressure or 
free-stream-disturbance amplitude of order E 2 h, i.e., and an induced-vortex spanwise 
velocity w of the size Eh 2 ln{h), i.e., E^ins, at most. 

The flow properties in the four decks are examined below in this section, for the main 
deck, in Section 3, for the buffer and lower decks, and in Section 4 for the upper deck. 
Since nonparallel-flow effects are incorporated for completeness, certain extra features arise 

including the basic-flow expansion profiles Xu b (y) ,v b (y) (i.e., u = u + (x — x 0 )u b H , v — 

£■;, + ••■) satisfying the steady boundary-layer balance u b + tJj = 0, u u b + v b t? = u", and 
A 6 = ul(0). 

In the main deck the unsteady 3D flow solution expands in the pattern 

u = fr(y) + £*Xu b {y) + eV 3) + erifZ* 1 ) + • • • + + e^llW + • • • , (2.2a) 

V = eT yf 1 ) + g: % + E^w + eM 3) + • • • + + • • • , (2.26) 

w = + e^W + Em ZE + e^u^ e) + • • • + e 5 uJ< 3) + • • • , (2.2c) 

p = H 1- e 6 ^ H (2. 2d) 

essentially. Here y = £*y, the superscript (3) denotes the induced vortex flow, which is 
independent of the fastest streamwise and temporal scales X,T, the superscript (1) refers 
to the main TS contributions present, while the double overbar denotes higher-order TS 
effects. The contributions and similar higher-order terms in u, v are forced by the 

induced-vortex flow outside the boundary-layer as explained in Section 4 below. Also, (e) 
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signifies the extra TS contributions affected by the nonparallelism (oc u b ,v b in (2.2a,b)). 
The governing equations for the successive terms above are obtained from substitution into 
the Navier-Stokes equations and are as follows. 

First, the dominant TS waves satisfy continuity and momentum equations which are 
quasi-steady and show the pressure variation across the boundary layer to be negligible: 


0$ + = 0, u = 0, (2.3a, 6) 

0 = -P^ , u fl#* = ~ P? • (2.3c, d) 

These yield the expected displacement-effect solutions 

a* 1 ) = {Ai\E\ + AwE^T? + c.c., (2.4a) 

= — ia.{A\\E\ -t- AwE^u -(- c.c., (2.46) 

P^ — (Pll-^l + Pu^) + c.c., (2.4c) 

= {wuEi + W 12 E 2 ) + c.c. : uTn i i 2 = z ff3p 1112 /2au, (2.4d) 


where c.c. denotes the complex conjugate and all the fast dependence on X, T is contained 
in the two oblique waves 


Ei, 2 = exp[t(aX ± fiZ/2 - flT)]. (2.5) 

Here the TS wavenumbers a,0 / 2 and frequency f 1 are real (see (3.14) below), and the 
displacement decrements An, ^12 and pressure components p n ,p 12 are unknown functions 
of X. 

Next, the higher-order TS influences represented by double-overbarred and (e) terms in 
(2.2) are controlled effectively by forced versions of the main TS balances (2.3) accounting 
for unsteadiness and nonparallelism. These yield in particular 

— u e \E\ -J- tT e 2 jE /2 4- c.c. (2.6 a) 


and so on, where 

u el = A el vf + AuXu'i, (2.66) 

v el = -(iaA e 1 + A n -^)u - iaA n Xu b , (2.6 c) 

with — A el representing an extra displacement correction which again is unknown in ad- 
vance. These influences are of little significance at this stage. 

The induced-vortex motion here in the main desk is also fairly passive, merely respond- 
ing to the efflux (oc j 4 33 ) from the buffer layer examined in the next section. Thus the 
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governing equations for u^ 3 ), t^ 3 ),ti> 3 £,p( 3 ) are analogous to those for the dominant TS wave 
but on the slower X scale and give the solutions 

u^ 3) = A 33 E 3 U 1 + c.c., = — A zz ^Ezu + c.c. (2.7a, 6 ) 

P^ = p 33 (X)i ?3 + C.C., U7 3 J5 = D$e(X)Ez/u c.c. (2.7c, d) 

where Z7 3 £ satisfies dD SE /dX = — i/?p 33 , the vortex pressure p 33 is an unknown function 
of X, as is the vortex displacement — /I 33 , and the spanwise dependence of the vortex is on 

E 3 = exp {i@Z), ( 2 . 8 ) 

Significant vortex-TS interaction does not occur in the main deck and is confined instead 
to the lower reaches of the boundary layer, specifically in the buffer and lower decks which 
are discussed in the subsequent section. The nonlinear interaction arises from inertial 
interplay which produces the combinations 

E\E^ X = Ez, EzEz = E\, E z x Ei = Ez, (2.9) 

thus enabling the two oblique TS waves to force the vortex motion, and the vortex, com- 
bined with one TS wave, to force the other TS wave, at the amplitude-squared level. In 
addition, the enhanced vortex-TS interaction near the wall is due to the increase in the TS 
spanwise motion there as represented by the singular behavior of w ^ in (2.4) as y — ► 0+, 
as well as by the higher-order TS effects, a singular response which is smoothed out in the 
viscous lower deck. 


3. THE BUFFER DECK AND THE LOWER DECK 


(a) The buffer. 

This extra, buffer, deck of width 0 ( 52 ) in y is brought into action because it controls 
the major response of the induced vortex flow, due in turn to the algebraic decay of the 
3DTS waves on exit from the lower deck, i.e. their growth in the lower extremes of the 
main deck. The buffer deck has the underlying pattern 


w 




11 5 15 A (A 

s * u H h e * vS ’ + • • • , 

(3.1a) 

■ • + H , 

(3.16) 

+ e~* i + • • • , 

(3.1c) 

) + • • • , 

(3. Id) 
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where now y = e?y, these expansions being implied mostly by the behavior of the solu- 
tions in the main deck for small y but partly also by the lower-deck features discussed 
subsequently. 

The governing equations for the dominant TS disturbance t^ 1 ), etc., are merely the 


continuations of those in Section 2, in fact, so that 

u^) = A (AiiEi + A 12 E 2 ) + c.c., = — iaXy^AnEi + A 12 E 2 ) + c.c., (3.2a, 6) 

pi 1 ) = ( pnEi + P 12 E 2 ) + c.c., u;! 1 ) = (3(—pnEx + P 12 E 2 ) /2Aay + c.c. (3.2c, d) 

with pn — p n , pn = Pi 2 - Likewise the nonparallel effect here is given simply by 

u (e) = \ A e + (XAj + 4 #) )^i, (3.3a) 

= — [A e x + A^Xy — AixXX^y — Aixu^ 3 \ (3.36) 


and the higher-order TS contributions continue the trend of algebraic growth, in the form 


u — /3 2 {PuEi + PnEi) / 4Xo?y + c.c., (3.4a) 

v = t'(ni4n - P 2 p n / 4aA - ap n /X)Ei + *(nA x2 - P 2 pu/4ctX - ocpi 2 /X)E 2 + c.c., (3.46) 

w — fith^/aAy. (3.4c) 

These last contributions play a vital role nonlinearly, however, in the generation and evolu- 
tion of the longitudinal vortex motion, through the major inertial response which is given 
by 

u^tbx + + v^wy + 

in the spanwise momentum balance, i.e., by Ezpnp^PiAot 2 — /? 2 )/4a 2 A 2 y 2 , from (3.2), 
(3.4) along with (2.9), as far as the vortex is concerned. The lower-order contribution 
cancels out, we note, and the logarithmic vortex effect w^\ etc., is only 
a passive one, discussion of which can be deferred for now. Thus the vortex motion of 
concern is controlled by the forced viscous-inviscid equations 


&33X + ®33y + *7?W 33 — 0, 


(3.5a) 


Xyu 33r + At> ss = u 33 yy, (3.56) 

X y™ 33 X + - /? 2 /40£ 2 )phP^ 2 / A 2 y 2 = w 3m , (3.5c) 

where we split u l 3 ) — ► u 33 E 3 + c.c. and so on. The boundary conditions on this vortex flow 
are to match, respectively, with the main-deck behavior for large y and with the lower 
deck for small y: 


u 33 — *■ A.i4 33 , v 33 A>i 33 j-y, tb 33 ay 3 , as y — ► 00 


(3.5 d) 
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u 33 = v 33 = 0 at y = 0, tu 33 ip( 1 - /3 2 /4a 2 )\ 2 p ll p > n iny + $(X) as y — >■ 0 + . (3.5e) 

Here the term $ in w 33 in (3.5e) is unknown as yet and has to be determined by the lower- 
deck response, as examined below, while the logarithmic term forces corresponding effects 
c x yiny, y z tny in t> 33 , u 33 as y tends to zero. All these effects merge with those of the 
lower deck. The main feature that the induced vortex motion imposes on the lower-deck 

a 

behavior is the skin-friction modification, u 33 ~ A 33 y as y — ► 0+ where 


A 3 3(X)e=<9 & «3 3 (*,0), (3.6) 

but that depends on the TS-induced forcing PnPj 2 (3.5c) and $ in (3.5e) which in turn 
depend on the skin-friction modification A 33 (see below), so that the whole process is highly 
interactive. 

Further points of note here concern: the extra logarithmic effects in (3.1); the fact that 
viscosity affects the vortex motion in the buffer zone; the pure displacement effect in (3.5d) 
and the feature that the pressure feedback from the upper deck due to A 33 E 3 (+ c.c.) is 
negligible in (3.5b, c); and how the entire structure might be affected by pressure gradients 
or cross-flow. These points are addressed subsequently. 

The logarithmically larger vortex effects and corresponding terms in u, v satisfy 
(3.5a-c) again, but with 33 replaced by 33L when the split tb$ — ► w 33 lE 3 + c.c. is applied 
and with the forcing term oc PuPi 2 absent. The appropriate boundary conditions then 
are w 33L = 7 ,u 3 3L = v 33L = du 33L /dy = 0 at y = 0 and w 33L -»• 0, du 33L /dy -*• 0 as 
y — ► 00 , where the sole forcing term 7 = i(3(l — / d 2 /4a 2 )A _ 2 p 11 pJ 2 /2 because of the lower- 
deck spanwise response in (3.18d) below (with Y — ► e~5y) and the zero-shear condition 
on u 33 l is required to merge the buffer- and lower-deck induced shears. The solution here 
is de-coupled in effect from the dominant TS-vortex interaction. Thus the solution for 
w 33L can be written formally in terms of a Fourier or Laplace transform in X and an Airy 
function, giving the required decay at large y, while the shear t 33 l = du 33 L/dy is governed 
by the same equation as is w 33 i except for a forcing \i(3w 33 L , and subject to f 33 i vanishing 
at y = 0 and as y — > 00 , the transform solution for f 33L therefore being a combination of the 
Airy function and its derivative. Integration with respect to y, subject to u 33 l vanishing 
at y = 0, then yields u 33 l for all y and the displacement correction, A 33L say, emerges from 
the large-y value of u 33 l- Hence the logarithmic (£ne) parts of the induced vortex motion 
here play only a passive part. 

b. The lower deck 

The lower deck, of width 0(e: 5 ) and containing the critical layer, lies beneath the buffer 
and brings in viscous forces to reduce the TS flow to rest at the wall. This smoothing of 
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the TS distributions in turn serves to smooth out the forced vortex response, in particular 
modifying the logarithmic irregularities connected with (3.5e) for small y. In the lower 
deck, where y = e 5 Y say, the flow solution can be expressed in the form 


u — eXY + e 2 [A 3 y + -|- s 2 «<*) + e 4 •+• 0(e*y 2 ) + • • • , (3.7a) 

v = + e — e 6 A 4 y 2 /2 4 , (3.7 b) 

xv — e^w^ + e^w^ + + • • • , (3.7c) 

p = 4 1- e 6 p( 3 ) + • • • , (3.7d) 


with Ag denoting X 33 E 3 + c.c. The induced-vortex components are principally the shear 
contribution X 3 Y (produced via the buffer) and the more complicated parts u( 3 ), u( 3 ), u>( 3 ), 
whereas the main TS components have superscript (1) again. The continuity, streamwise 
momentum and spanwise momentum balances then yield, successively, and in view of the 
interplay in (2.9), 


< ag* + v? + 4 1} 

( »!?' 1 + 4 3) = 0; 


0, uty + 4- vfi + 4 e) = 0, 


(3.8a — c) 


f + AVu^ + t/^A = — Pj^ + ttyjr, 


4 e) + Ay(u^ + ug. } ) + (XA 6 y + A 8 y)ttg ) + sWfA* + Xa 4 ) + »«a 
- 4 e) aii) + r.w 

[ + u( 3 )A 4- 4- u^ 1 ) 4 1} = 


(3.9a - c) 


4 1} + Ay 4 1} = -pg* + 4 1 y> 

■ 4 e) + Ay (4 } + + (^Ajy + A 3 y)wg } = -pg } + 4r, (3.10a - c) 

uWwg) + = Wyy, 

for the main TS, the forced TS and the vortex flow, in turn, all three of which need to be 
addressed below. Only the E 3 components are of concern in (3.9c), (3.10c), we note, and 
the y-momentum balance requires the pressure effects shown in (3.7d) to be independent 
of y, leaving the results p^ 1 ^ = pW = p^, p^ = p M = pW,p( 3 ) = p( 3 ) = p^ 3 ). 

Splitting the TS responses into their pairs of oblique waves and the vortex into its 
spanwise parts now, so that 


u (1) = u n Ei + u 12 I?2 + c.c., u {e) = u e \Ei + u e 2 E 2 + c.c., u (3) = U 33 E 3 + c.c. (3.11) 
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and so on, we have from (3.8a), (3.9a), (3.10a) the governing equations 


taU ln + Vi n Y ± ifiwin/2 = 0, 

-iUuin + ia\Yu ln + Avi„ = -iap ln + u ln yy, 
— in Win + iaXY U>in = ^i/3p ln /2 + WinYY, 


(3.12a) 

(3.126) 

(3.12c) 


for the main TS wave, for n — 1,2 respectively, subject to the constraints of no slip at the 
wall, of matching with the buffer solution and of pressure-displacement interaction: 


«ln = Vin = W\ n = 0 at Y = 0, (3.12 d) 

Ui n — *• A-A ln , win —*■ 0(Y _1 ) as Y —* oo, (3.12e) 

(a 2 + 0 1 / 4) bm = « 2 ^m- (3.12/) 

Here (3.12f) anticipates the results of Section 4. The solution of the linear TS problem 
(3.12a-f) can be expressed in terms of the Airy function At and leads to the eigenrelation 

( Co = 

\ 2 Ai'(£o)/K = (taA)^(a 2 + (3 2 /4)*, with < (3.13) 

l K = J£Ai( s ) ds > 

determining a for given real /?, fl values. The required neutral condition corresponds to 
£o = — dit’r, Ai'(£o)/ K = dii^\d\ « 2.30, d^ « 1.00], however, so that 

a® (a 2 + /? 2 /4)® = d^X*, fl = diCc^A® (3.14) 

relate the real a,/?, fl values of concern here. In the above i* denotes exp(tV/6). 

Next, the forced TS wave (u e i, etc.), driven by the combination of modulation and 
nonparallel- and vortex-flow effects in (3.8b), (3.9b), (3.10b), satisfies (3.12a-c) again with 
subscripts 11 replaced by el of course and with u U jf, XYu n x + XX b Y iau u + A 33 Y to:ui 2 + 
XCXbVn + A 33 vi 2 + Pnx’ AYw llT + XX b Yiawn + X^Yicnwn added to the respective left- 
hand sides, for n = 1, and similarly for n = 2. In addition, the boundary conditions are 
(3.12d,e) again with 11 — *• el and with X b JCA n + A 33 A 12 supplementing the right-hand side 
of (3.12e) (for n = 1, similarly for n = 2), coupled with the interaction law 

(a 2 + f3 2 /4)b en = ~ 2 taAjjp + ia 3 A in x/ (a 2 + /? 2 /4) (3.15) 

stemming from the upper-deck properties presented in the next section. This forced- 
TS system requires a compatibility relation to hold between all the driving effects just 
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described, and, after some manipulation, the relation is found to provide the governing 
equations 

f a Piix + bXXiPn + c\ 33 X- 1 p n = 0 , 


{ a,Pux + b XX b p 12 + cA; 3 A 1 p n = 0 


(3.16a, 6 ) 


for the slower-scale modulations of the primary TS pressure amplitudes. In each of 
(3.16a,b) the second term represents the effects of the nonparallelism of the basic flow 
while the third term comes from the extra induced-vortex motion (via (3.6)). Also, 


a = 2Bt 0 riD/3a&? - iXB~*(4/3 + 0 2 /12a 2 ), (3.17a) 

b = -2ai 0 r x BDl3^ - 5B*/3a, (3.176) 

c = iDZoA-*{2B ri /3 + p 2 r 2 / 8 ) - XoT 1 B-*{$cl 2 /3 + /? 2 /24), (3.17c) 

where D = l + /c£ 0 /A*'(£ 0 ),.B = a 2 + /? 2 /4, A = iaX,r x = At(£o)M*'(^o), r 2 = /c/At(6)- An 
origin shift in terms of X is implicit in the above to account for the sub- or super-criticality 
of the flow. 

Finally here, the induced vortex flow in the lower deck is controlled by (3.8c), ( 3 . 9 c), 
(3.10c) with (3.11). But the spanwise balance (3.10c) alone is sufficient in fact to pin down 
the most significant unknown property, the viscous TS-amplitude-squared contribution 
(or jump term) $(X) which is required to complete the prescription of the induced vortex 
motion in (3.5). Thus, taking the E 3 component in (3.10c), we have to solve 


w'ss = + iau* 12 w n + v n w 22 + v* 12 w' n + tpw* 12 


(3.18a) 


for {^(X, Y”), subject to no slip at the wall and to joining with the buffer deck: 


u >33 = 0 at Y = 0 , (3.186) 

1633 — ► 0 as Y — ► 00 , (3.18c) 

=» u >33 ~ -*/?(! - p 2 /4a 2 )X- 2 ? n p[ 2 tnY + $(X) as Y -> 00 . (3.18d) 

The logarithmic behavior here, arising from the algebraic decay of the 3DTS velocities 
involved on the right of (3.18a) into their larg e-Y asymptotes (from (3.12)), matches with 
the corresponding logarithmic responses present in (3.1), (3.5), and similarly for u S3 ,v 23 . 
Further, in view of the nonlinear forcing in (3.18a), the jump term $ can be written 

* = # 11^2 (3-19) 

where <f> is a constant to be determined, for prescribed /3. Relevant comments at this stage 
are the following. First, we need the complete behavior 0 (tnY) + 0 ( 1 ) in (3.18d) to feed as 
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a boundary condition into the main vortex problem (3.5) in the buffer. Thus, the whole 
balance between the buffer and lower decks is interactive. Second, some of the solution 
here follows by taking a limit of the Hall and Smith (1984) analysis to deal with standing 
waves. Third, a nonzero basic pressure gradient adds an e 2 Y 2 term to u, which is a minor 
effect on the vortex; see also Section 6. Fourth, the pressure e 6 ps is induced through the 
upper-deck interaction, due to displacement: see the next section. Fifth, the effects of the 
induced vortex and of the basic-flow nonparallelism are comparable here and in the buffer. 


4. THE PRESSURE-DISPLACEMENT INTERACTIONS VIA THE UPPER 
DECK 


The upper deck lies outside the boundary layer and has a y-scale of 0(e 3 ), within 
which linearized potential-flow properties act to link together the various pressure and 
displacement distributions present. The velocity and pressure fields here have 


(u,v, w,p) = (1,0, 0,0) + e‘* (u^, + e «*(u^, v^ e \w^ e \p^) 

+ (e 6 u^ 3 \eJv^ 3 \e^it;^ 3 \e 6 p^^) + • • • , (4.1) 

where the terms on the right-hand side represent, in order, the uniform external stream, 
the main TS wave, the extra forced TS wave and the induced vortex. The disparity in 
the vortex velocity sizes is due to the relatively slow streamwise variation of the vortex 
motion. Also here, y = e z y' with y' of order unity. 

So the main TS wave is governed by the linear inviscid system 


ttj? + + w z ] = 0> = -Pjp>Vjp = -P { y\ wjp = ~Pz\ (4-2) 

p^ — ► p^, — ► ot z (AuEi + A 12 E 2 ) + c.c. as y' — ♦ 0+ (4.3a, 6) 


and with suitable boundedness conditions in the farfield. The constraints (4.3a, b) match 
the pressure and normal velocity with those in the main deck of Section 2, thus ensuring 
the matching of the streamwise and spanwise velocities as well. A decomposition as in 
(2.4) applied to p^ then produces the pressure-displacement law anticipated in (3.12f), 
since pM satisfies Laplace’s equation, from (4.2). 

In similar vein, the forced TS contribution satisfies in effect a forced version of the sys- 
tem (4.2) - (4.3) and that leads to the law quoted already in (3.15). The induced vortex flow , 
on the other hand, is governed by the effective cross-flow balances 


«L 3) + 4 3) = 0 , 


(»)__«(») ...m -_„(?> 


= ~P 


r' J 


W 


x - p z 


(4.4) 
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here because of its comparatively small u,p parts, and the streamwise balance u^- = —p^- 
is subsidiary. Hence under decomposition, and merging with the main-deck pressure, 

P (3) = P33 exp(-/V)-Es + c.c. (4.5) 

and the displacement-or v-match then gives the result 

PPzz — ~A-z 3X x- (4*6) 

This fixes the vortex pressure p 33 (X) in the boundary layer in terms of the displacement 
— .A33 which is to be determined by the buffer- vortex problem (3.5). We observe that the 
induced pressure e 6 pW above has relatively little feedback effect on the properties holding 
inside the boundary layer. Thus, the induced vortex’s effects on the TS waves, during 
the present regime, are confined predominantly within the boundary layer (and, further, 
within the buffer and lower decks) , even though the actual vortex motion produced spreads 
significantly far outside the boundary layer into the upper deck. See also Fig. 1. 

The upper-deck vortex pattern also produces nonzero slip-velocity components u( 3 ), ti/ 3 ) 
at the upper edge of the main deck, as y' — ► 0+. These induced slip velocities tie in with 
the expressions proposed previously for the main, buffer and lower decks, e.g., see the 
term tZ > 3E in (2.2), but, like the induced vortex pressure above, they play no role in the 
dominant nonlinear interaction between the TS wave and its induced vortex. This TS- 
vortex nonlinear interaction is addressed further in the following section. 


5. THE NONLINEAR INTERACTION EQUATIONS, AND SOLUTION 
PROPERTIES 


In summary, the nonlinear induced vortex/TS interaction requires us to solve the fol- 
lowing nonlinear partial-ordinary differential system consisting of three parts: first, 

r U Y +V t + ifiW = 0, 

< yU-jf+V = Uff, (5.1a — c) 

. yW Y +iK\P\ 2 y-* = W n , 


( U -► 0, V -*■ 0, W ~ —iK\P\ 2 tny + $(X), as y 0, 

[ U — * A, V Ay, W oc y -3 , as y — *■ 00 

for the vortex (17, V, W)\ second, 

P~x + (ciAjX + C2A33 )P — 0 


(5. Id, e) 


(5.2a) 
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for the TS amplitude (P), where 


X 33 {Y) = U i {X,0); (5.26) 

and, third, using the expression (3.18a), 

f Y { C 2 *33 (X,Y 1 )dY 1 }dY 2 ~ -iK\P\HnY + $(X), (5.3a) 

Jo j oo 

=► $(X) = |P|V (5-36) 

for the determination of the constant <£, and hence the function $, for a given spanwise 
wavenumber f3. Here U,V,W,P,A can be related back to U33, 033, *33, Pu (or p 12 ), A33, 
although to fix matters we are focusing attention now on the case where p n = p 12 and A33 
is real. Also, 

K = P{ 1 - (3>/4a 2 ), (5.4) 

while the constants C\ = 6/a, C2 = c/Aa in (5.2a) follow from evaluation of a, 6, c in 
(3.17) for given /? , and likewise for a from (3.14). The solution properties are addressed 
below, although most often with the U — V — W formulation in (5.1) replaced by a t — W 
formulation where W is still to be derived from (5.1c-e) but 

T vi ~ V t Y = (5.5a) 

t — ► 0 as y — »• 00, Ty = 0 at y = 0 (5.56) 

is then used to find r(= Uy) and hence 

A 33 (X) = t(X,0). (5.5c) 

Here (5.5a, b) follow from (5.1a,b,d,e) and U,V,t are real whereas W is pure imaginary. 

Computational solutions of the interactive problem in question, (5.1c-e) with (5.2a), 
(5.3), (5.5), were sought with a simple implicit finite-difference procedure marching forward 
in X, for several spanwise wavenumbers (3 defining K in (5.4) , as follows. Given the solution 
up to the station X = X n _i say, the W profile at the next station X n = X n _i + AX is 
obtained by inversion of the tridiagonal system 

(wtt+t -2K+wJj-i)/±y 2 -MK- w tij)/±x = -iiiiMin), i = 2 -» j-i, ( 5.6a) 

W nl = <t>\Pn-l\\ (5-66) 

W* = + lo in{yj) (5.6 c) 

representing (5.1c-e), after the substitution W = —i 0 tny + W A to keep W A finite at y = 0. 

Here AX, Ay are suitably small grid sizes, the subscripts n,j denote the function values 
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at the gridpoints X n , y,[= (j — l)Ay for j = 1 to J],7o = +tif|P| 2 ,yj is a suitably large 
edge value of y, and 7o,7i[= d^/dX] in (5.6a, c) are evaluated from the known previous 
solution values, as with |P n -i| 2 in (5.6b). Also, the number <f> is prescribed for a given 
/?: see below. Next, the r profile at X n follows in a similar manner from inverting the 
tridiagonal form 

(Tnj+ 1 - 2 T nj + T n j-i)/ Ay 2 - yj{r n j - V„)/Al = -i0{W$ - lotn{yj)} ( j 

T nl = T n2 , 

TnJ = 0 , 

associated with (5.5a, b), thus yielding the skin-friction factor (5.5c) as 

(A 33 )n = Tnl. (5.6ff) 

Finally at this station, P n is found from an analogue of (5.2a), 

P n = P n - 1 + AX[cxAjX n + c 2 (A M ) n ]JVi, (5.6/i) 

and then the scheme (5.6a-h) can be moved on to the next station downstream, and so on. 
Many other numerical procedures suggest themselves of course, but the above procedure 
proved to be stable and accurate for sufficiently small step sizes: see also Figures 2-5 where 
sample interaction results are presented, all starting at position X = — 1 upstream of the 
neutral TS point. 

The results shown are for the specific cases /? = 1,1.6, 2,2.5 [see reasons later], for 
which the constants involved in the nonlinear interaction equations are as follows (see also 
Figure 1(b)): 


(5.6d) 

(5.6e) 

(5.6/) 


p 

a 

0 

a 

b 

c 

4> 

1 

0.904 

28.85° 

-1.09 - 2.84i 

-0.59 - 0.99i 

-0.56 - 0.49i 

12.32i 

1.6 

0.754 

46.70° 

-1.59 - 3.66i 

-0.83 - 1.20i 

-0.70 + 0.17i 

63. 6i 

2 

0.617 

58.32° 

-2.53 - 5.24i 

-1.09 - 1.57i 

-0.84 + 0.88i 

179.3i 

2.5 

0.433 

70.89° 

-5.79 - 10.78i 

-1.75 - 2.51i 

-1.25 + 2.32i 

679i 


The normalized case A = 1,A* = — 1 is assumed here, while the a,b,c values follow from 
numerical evaluation of (3.17), a from (3.14) and <f> comes from (3.18) - (3.19) (or (5.3)) 
after some extensive grid-effect studies to allow for the rather slow approach to the asymp- 
tote there. We observe that since A 33 is real P could be replaced by |P| in (5.2a) with 
Ci r , c 2r then replacing Cj,c 2 respectively. Also in the table above 0 gives the TS wave angles 
tan _1 (/?/2a). 
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The ultimate behavior of the nonlinear interactive flow, as X increases, is considered 
analytically now. There appear to be several possible options, the first of which has an 
algebraic singularity arising at a finite position, say as X — ► X 0 —. The orders of magnitude 
suggest the balances r ~ A33 ~ |djjf|, from (5.2a), (5.5c), |d & | 3 ~ \d^\ and W ~ \dx\* from 
(5.5a), and hence P ~ approximately, from (5.1c,d), (5.3). So, with allowance for 

the logarithmic term in (5. Id), the proposed local response as X — ♦ X 0 — is of the form 

|P|~(X 0 -X)-*ir*|P|, W ~(X 0 -T)-iw, (5.7a, 6) 

(r, A33) - (X-X) _1 (r,A33), r, = y(X 0 -X)-£, (5.7c, d) 

to leading order, with L = —ln{X 0 — X) and P, A 33 unknown constants. Here the vortex 
system reduces to the ordinary differential equations for fV(tj), r(tj), 



W" - rfW'/Z - 5r]W/Z = 0 

(5.8a) 

subject to 

f" - rj 2 ?/ 3 — r)T = -ifiW 

(5.86) 

W{< 

x>) = f(oo) = r'(0) = 0, W’(O) = tP:|P| 2 /3 

(5.8c - /) 

and the TS part requires, 

for consistency, 

5 - 

(5.9) 


- = — C 2 rA33 
6 


where A33 = f(0). It is worth noting that here and in the subsequent option the TS forcing 
on the vortex makes itself felt only through a constant inner constraint on W , here (5.8f), 
which is due to the logarithmic contribution in (5. Id), the other forcings [oc K in (5.1c) 
and <f> in (5. Id)] being negligible at this level, as is the linear-growth term oc c\ in this case. 
The solution for W as a function of a = rj 3 is 

W(a) = K , /” r»(r + i yh'^dr, K, = .if |>|V{3*r(|)r(|)} (5.10a) 

and therefore 

t(o) = o*t a l q f 3 e~ ,l/q {-ip f M^(«2)s2 *ds 2 + B}dsi/9 (5.106) 

J oo J oo 

with B chosen to make dr /dr) zero at a = 0, from (5.8a-f). This then gives A 33 = 

„ _ I 

— [i(3 iy(s 2 )s 2 3 ds 2 + P]/3. Hence the existence of this local algebraic break-up depends 
on the relation 

c 2r K\P\ 2 = d 3 (5.11) 
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being satisfied, because of (5.9). Here d 3 is a positive real constant. Since |P| 2 is to be 
positive, the break-up (5.7) is therefore a candidate for the ultimate interaction behavior 
in the regime where c 2r K is positive. This regime covers all wave angles 6 except for a 
small interval of 0 between 45° and approximately 50° where Kc 2r goes negative, as shown 
in Fig. 1(b). 

The second option, in contrast, has the solution continuing to downstream infinity, 
where an algebraic response takes effect, namely 


|P| ~ X«L"5|P|, W ~ X*W{r}), 

(5.12a, b) 

(r, A 33 ) ~ X(r,A 33 ), rj = 0C*, 

(5.12c, d) 

as X — * oo, with L = inX now. This option arises from balancing the growth terms 
oc c l5 c 2 in (5.2a), yielding (5.12c), then (5.12b,d) follow from (5.5a), after which (5. Id) 
suggests the size in (5.12a). Thus here the X-derivative in (5.2a) and the forcing terms 
oc K, $ in (5.1c,d) are negligible to leading order. The controlling equations and conditions 

for this option are 


W" + T] 2 W'/3 - nW/3 = 0, 

(5.13a) 

t" + r) 2 P/3 — r)T = —i0W, 

(5.136) 

W(oo) = r(oo) = f'(0) = 0, 1^(0) = —iK\P\ 2 /3. 

(5.13c) 

0 = c lr A 6 + c 2r A 33 , where A 33 = r(o), 

(5.13d) 

from substitution of (5.12) into (5.1)-(5.5), and so we have the solutions 


W = i(K\P\ i /3)r] [* rji 2 exp(— I7i/9)dr7i, 

J OO 

(5.14a) 

^ poo r oo 4 

f = P(K\P\ 2 /27) / s i 3 exp(— si/9)dsids 2 . 

J a J $2 

(5.146) 


Here a = r/ 3 . Therefore A 33 = /?(.ff|P| 2 /9)r(2/3)33, from (5.14b), and (5.13d) then gives 
the relation 

0\P\ 2 c 2r K = d 4 (5.15) 


where, since A* = —1 and c ir is always positive, the constant di is positive. In consequence 
this option 2 also applies for positive values of Kc 2r only. 

Option 3 is that the nonlinear solution amplifies exponentially as X — > oo, with a 
growth exp(sX/2), say, in |P| forcing a vortex flow of amplitude exp(sX) (from 5.1c,d)) 
which in turn requires the c 2 term in (5.2a) to dominate alone. Thus 

|P| ~ exp(sX/2)|P|, W ~ exp(sX)jy, (5.16a, 6) 
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(5.16c, d) 


(r,A 33 ) ~ exp(sX)(f, A 33 ), y ~ 1. 

Then from (5.1)-(5.5) the leading balances are 

f W" - syW = iK\P\ 2 y~ 2 , 

< (5.17a, 6) 

t" - syr = -i/3W 

with 

l^(oo) = f(oo) = r'(0) = 0, W ~ \P\ 2 {—iK iny + <t>} as y 0 (5.17c - /) 


and 

A 33 = ?(0) = 0. (5.18) 

The positive constant s is to be determined. The solution of (5.17) can be expressed as 

W = iK\P\ 2 M(y), r = s~i(3K\P\ 2 N{y), y = s*y (5.19a) 

with 

M{y) - BiAi(y) - iny - Ai(y) £ Ai~ 2 {y{]{j ^ y 2 Ai(y 2 )iny 2 dy 2 }dy l , (5.196) 

N(y) - B 2 Ai(y) + M'(y) + y~ l + Ai(y ) J Ai~ 2 (yi){ J Ai(y 2 )dy 2 }dy l (5.19c) 

where Ai is the Airy function again and BiAi(0) = (tns)/3 — i<f>/K,B 2 Ai'(0) = l/(3At'(0)). 
Hence A 33 oc 7V(0) follows, and the requirement (5.18) leads to the real positive value 

s = exp(— Ci/ e 2 ) (5.20) 

for the growth factor s, where ei = 1.21 — 0.73t<f>/K, e 2 = 0.243. This option, in which 
(unlike in options 1,2) all the TS forcings on the vortex flow play a substantial part, is 
therefore a versatile one in the sense that it is available for all values of K, c lr , c 2r , i.e., for 
all wave angles. 

Option 4 has the TS disturbance |P| becoming very small/negligible, and the vortex 
flow then grows slowly on its own with distance X, from its initial state upstream. This 
option is ultimately unstable to the TS waves (if present) however, via (5.2a), since the 
vortex skin friction A 33 is insufficient to counteract the nonparallel-growth term oc Aj 
downstream. 

The computations described earlier tend to support option 1 in the cases where Kc 2r 
is positive, i.e. for the spanwise wavenumbers /3 = 1,2, 2.5, although the finite-distance 
break-up associated with that option can be delayed for a very long distance, depending 
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on the starting conditions upstream: see the figures. It is interesting that the long delay 
corresponds in fact to the near-attainment of option 4. For the case (3 = 1.6, in contrast, 
where Kc 2r is negative, the computations seem to point to the long-scale exponential 
option 3 being approached downstream, the rapid growth in that case being ascribed to 
the enormous nonlinear growth factor s « 10 410 7 from (5.20), which in turn is due mostly 
to the numerical largeness of \<f>/K\ when (3 = 1.6. The value of s is also large for the other 
spanwise wavenumbers although there the option 1 appears to take preference, as might 
be expected. The relevance of all these options (and we note that no others have been 
found yet) is discussed in the next section. 

6. FURTHER DISCUSSION 

Two immediate extensions of the theory should be mentioned first and they concern 
external pressure gradient effects and unequal TS amplitudes, respectively. The influence 
of an external pressure gradient or non-uniform slip velocity just outside the boundary 
layer is accommodated by the present nonparallel-flow theory through the term involving 
ci in (5.2a), which in turn comes from the normalized TS dispersion relation (3.13) in effect. 
In dimensional variables (denoted by “D”) this shows that the frequency Qd is related to 

3 1 _ 

the triple-deck frequency, Qtd say, by f l# = (duD/dyD)lfiTDU 0 oi / oo{'U'D)7 2 locally, where 
w stands for the dimensional value at the wall and oo for the dimensional reference value, 
e.g., far upstream, while (u£>) e is the external slip velocity. For a fixed input frequency 
f Id upstream (and fixed Uoo, !/<„), therefore, f Itd increases or decreases, i.e. the TS 
waves are destabilized or stabilized (corresponding to A j = =Fl i n effect), according as the 
combination (dup / dt/r>)u>(up)7 2 decreases or increases, in turn. Thus the issue depends 
on the basic boundary-layer slip velocity and the wall shear. For instance, a favorable 
pressure gradient has (uu)„ increasing but that then tends to increase the wall shear also, 
so reducing the stabilizing trend. A sufficiently favorable presure gradient could lead to an 
effective change in the sign of X b , nevertheless, in (5.2a), and that then alters the nonlinear 
interaction properties, as indeed the options 1-4 in Section 5 would suggest. An example of 
a computation associated with the stabilizing value A t = 1 is shown in Figure 6. Next, the 
special case of equal TS pressures p n , p 12 , taken in Section 5, holds also for the broader case 
where |p n |, |p 12 | are equal with the phases of p n ,p 12 differing, i.e. the phase difference is 
insignificant then. The most general case of unequal |p u |, |p 12 | however is not covered yet. 
It might be argued that the options 1-4 describe the ultimate response even in the most 
general case but the effects of unequal amplitudes and the associated rise in significance of 
the phases remain to be seen. The special case of 45° waves where the interaction constant 
K is zero may require more research as well. 
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Comparisons between the present study and corresponding experiments as well as pre- 
vious computational work tend to be favorable on the qualitative features, if not more. 
This is particularly so with respect to the implied abrupt break-up or change in scale of 
the streamwise vortex motion (in option 1), due to the generation of the oblique TS waves 
and then nonlinear interaction. A number of other features (e.g., see scales below) also 
seem common to the theory and the experiments. 

The present approach is based on a rational scheme throughout, we should stress, which 
is aimed among other things at shedding light on the scales which are central to TS/ vortex 
interaction and hence to transition. The scales of the pressure and velocity disturbance 
levels associated with the current nonlinear-interaction regime for instance are quite small, 
as (2.2), (3.1), (3.7), (4.1) show; e.g., the streamwise velocity disturbance implied in (4.1) 
is approximately 0.4% of the free stream velocity at a global Reynolds number of 10 6 . The 
theory also accounts for nonparallel-flow effects in a rational way, which is a vital aspect 
in the boundary-layer context, and as yet there appears to be no alternative acceptable 
theoretical approach. 

Addressing the further implications of the theory, we feel that the main repercussions 
stem from the options 1-4 for the terminal behavior of the nonlinear TS/induced-vortex 
interaction described at the end of the previous section. See the sketch in Figure 7. 
Option 1 indicates a shortening of the streamwise interaction scale, as the break-up in 

(5.7) takes place, together with enhanced vortex and TS amplitudes. In fact, when Xo — X 
shortens to 0(et) the reduced pressure amplitude \P\ rises to 0(e - i), formally, apart from 
a logarithmic factor, thus forcing the unsealed wall pressure p (e.g., in (3.7)) to rise to 
0(e 2 ). Along with that, the scaled spanwise vortex velocity W increases to 0(e“l), from 

(5.7) , meaning that the un-scaled vortex w-contribution increases to 0(£*.£ - *), i.e. 0(e), 
from (3.1) or (3.7), in the buffer or lower deck; and the typical buffer thickness 0(e!|y|) 
becomes O(et.e^), i.e., 0(e 6 ), coincident with the lower-deck thickness. All these scales, 
and the others present, imply that the next stage encountered beyond the break-up in 

(5.7) brings in the fully nonlinear unsteady 3D triple-deck response on the 0(e 3 ) length 
scale (streamwise and spanwise), that is, the 3D unsteady boundary-layer equations with 
pressure-displacement interaction. The 3D triple-deck description (see computations by 
Smith (1988)) then controls both the nonlinear TS and the induced-vortex behavior. An 
interpretation of the above is that the nonlinear 3D TS/vortex interaction addressed in 
the present work provides a powerful mechanism for small input disturbances upstream to 
“burst” abruptly through to full nonlinear-TS status, at a sub- or super-critical location, 
via the singular behavior in (5.7). 

Option 2, on the other hand, has the streamwise scale lengthening downstream as 
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indicated in (5.12). Indeed, when the typical un-scaled distance x enlarges to the airfoil 
scale 0(1), i.e. X — > 0 (e: - 2 ) formally, then the vortex skin-friction factor A 33 rises to 0(e“ 2 ), 
from (5.12), suggesting an increased contribution of order e in the lower-deck velocity u 
in (3.7a). This is comparable with the basic-flow contribution (erAT) there and so changes 
the TS stability characteristics. Simultaneously the buffer thickness 0(e:t|y|) rises to 0(e 4 ), 
meaning that the main boundary layer is affected more substantially than before by the 
vortex flow, and the u-velocity component of the vortex rises to 0(1). Thus the entire 
boundary-layer flow is also altered, in the implied next stage. This long-scale next stage 
is examined by Hall and Smith (1988). 

The subsequent or “second” stages inferred above, as well as the long-scale ones cor- 
responding to options 3,4, clearly merit further study. This is partly to understand more 
the ensuing behavior in the TS/ vortex interactions which start in the present regime (of 
(2.1)), of course, but partly also to tackle the interactions that can arise with higher input 
amplitudes leading directly into those second stages, i.e. the by-pass process. That process 
makes all the second stages available in effect, provided the latter are self-consistent as in 
options 1-4, even if they cannot be reached via the smaller-amplitude first stage. 

Another aspect worthy of further investigation corresponds to taking the spanwise 
wavenumber 0 small, to study the 3D de-stabilization of an incident 2D nonlinear TS 
wave. The scales in the present regime provide a starting point for such an investigation, 
including the influence of nonparallelism. The theoretical approach is believed to apply to 
a number of other interesting flow configurations as well, and to have a strong connection 
with some observed transition phenomena such as the creation of “Lambda vortices.” 

Finally, the importance of nonparallel-flow effects in certain interaction cases is re- 
emphasized. Some remarks on them have been made already in this section. Nonparal- 
lelism of the boundary layer is vital in both of the options 2,3, corresponding to long-scale 
induced vortex flow, and hence also in the second stages subsequent to those options. Ne- 
glect of the nonparallelism seems acceptable for the shortening interaction of option 1 and 
the like but it would lead to erroneous results otherwise. 
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FIGURE CAPTIONS 

Figure 1. (a) Diagram of the 3D TS/vortex structure for nonlinear interaction in a 2D 

boundary layer (BL): see also Fig. 7 below. 

(b) Certain coefficients plotted against the TS wave angle 9. 

Figure 2. Nonlinear-interaction computed results for 0 = 1, wave angle 9 = 28.95°. 

(a) Vortex shear A 33 vs. X, (b) TS amplitude squared (|P| Z ) vs. X. Initial 
|P| = 0.1, initial vortex [ W,t ] = [qiy,qt{l + y 2 )] exp(-y z ). Case I: (91,92) = 
(0,0). Case II: (91,92) = (1,1). Grid has (AX, J, Ay) = (0.002,801,0.08) in 
general. Results 0, X are for grids (0.001, 801, 0.08), (0.001, 801, 0.04) in 
turn. Similar trends resulted for starting values |P| = 1,0.01,0.001. 

Figure 3. Computed results for 0 = 1.6, wave angle 9 = 46.70°. Starting |P| = 0.1, 
with (91,92) = (— 5, — 5), (1, 1), (5,5) for cases I, II, III respectively. Grid 
(0.001, 1601, 0.08); doubling AX, Ay produced little change. Results for 
case (0,0) of starting vortex at rest are very close to those of case II. 

Figure 4. Results for 0 = 2, wave angle 0 = 58.32°. Start |P| = 0.1, (91,92) = (0,0); 
results for (1,1) are almost identical. Grid (0.001, 1601, 0.08). 

Figure 5. Results for 0 = 2.5, wave angle 9 = 70.89°. Start |P| = 0.1, (91,92) = (0,0); 

grid (0.002, 1601, 0.08). Results for starting |P| = 0.1(») and for (1,1) (X) 
are also shown, as is typical effect (0) of doubling/halving AX, J, Ay. Similar 
trends resulted for other starting conditions calculated. 

Figure 6. Results for stabilizing case X b = 1. 

Figure 7. Summary of (a) effect of wave angle of nonlinear breakdown and (b),(c) next 

higher-amplitude stages of short/long scale following breakdown options 1,2. 
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Fig. 1(b). 
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